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The quest of novel compounds with special structures and unusual functionalities continues to be a central 
challenge to modern materials science. Even though their exact structures have puzzled scientists for 
decades, superhard transition-metal borides (TMBs) have long been believed to exist only in simple crystal 
structures. Here, we report on a polytypic phenomenon in superhard WB 3 and MoB 3 with a series of 
energetically degenerate structures due to the random stacking of metal layers amongst the interlocking 
boron layers. Such polytypism can create a multiphase solid- solution compound with a large number of 
interfaces amongst different polytypes, and these interfaces will strongly hinder the interlayer sliding 
movement within each polytype, thereby further increase the hardness of this particular material. 
Furthermore, in contrast to the conventional knowledge that intrinsically strong chemical bonds in 
superhard materials should lead to high lattice thermal conductivity, the polytypic TMB 3 manifest 
anomalously low lattice thermal conductivity due to structural disorders and phonon folding. These 
findings promise to open a new avenue to searching for novel superhard materials with additional 
functionalities. 
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Superhard materials are technologically important in many applications, from reducing the wear of everyday 
objects to creating machining tools. The hardness of a material is usually measured by indenter techniques 
and is thus determined by how resistant the bonds in the material are against distortions and how disloca- 
tions are able to move in the system 1 . A three-dimensional strong covalent network formed by carbon atoms 
makes diamond the hardest known material 2 . To overcome the limitation of diamond to ineffectively cut ferrous 
metals, dense TMBs with high boron content (e.g., OsB 2 , ReB 2 , CrB 4 , FeB 4 , MnB 4 ) 3 " 9 have been pursued. Implicit 
in such efforts is the assumption that dense TM atoms (creating high valence-electron density) prevent the 
structures from being squeezed while the high boron content (forming strong covalent B-B and TM-B network) 
withstands both elastic and plastic deformations, both of which enhance the hardness of materials. Amongst these 
TMBs, the highest boride of tungsten has currently become a focus of promising superhard materials because it 
combines ambient pressure synthesis, inexpensive raw materials with a measured superhardness (43.3-46.2 GPa 
under an applied load of 0.49 N) 10 " 17 . However, until recently, its structure has not been completely resolved. The 
heavy tungsten atoms hamper the accuracy to which the light boron atoms can be located by x-ray diffractions. 
Moreover, the versatile nature of boron atoms to form sp-, sp 2 -, and sp 3 -hybridized bonds brings about the 
coexistence of various phases. These technological challenges result in the structural and compositional uncer- 
tainties. Combining first-principles calculations with thermodynamics, the long assumed WB 4 phase 18 has 
recently been determined as a simple and stoichiometric WB 3 structure (called 2H) 19 . This conclusion has 
subsequently been substantiated by many theoretical and experimental studies 20 " 26 . On the other hand, agreement 
is far from complete for WB 3 . A theoretical study of WB 3 has suggested that the 2H structure is incompatible with 
experimental results (e.g., superhardness, normalized cla ratio) and that the question of the crystal structure of 
samples must be reopened 27 . In particular, an energetically more favorable WB 3 (called 3R) has been proposed 28 . 
Hence, it is still unclear whether other phases are possible and what structural similarities may exist amongst these 
phases. 

Since molybdenum is isoelectronic with tungsten, the crystal structure ambiguity for molybdenum borides is 
very much in parallel to that for tungsten borides. The highest boride of molybdenum was first reported as the 
MoB 4 phase 29 and later identified to be a triboride-based Mo!_ x B 3 structure 30 , but it was recently proposed 31 that 



SCIENTIFIC REPORTS | 4:5063 | DOI: 1 0.1 038/srep05063 



1 



H 



(a) (b) 

Figure 1 | Structures of boron and metal layers, (a) Overlap view of hexagonal boron (H) and metal layers (A) along the [001] direction, (b) Side views of 
boron (H) and metal layers (A, B, C) along the [110] direction. The honeycomb metal layer (A) can be derived from a close-packed metal layer by 
removing one third of the metal atoms (marked by the red plus sign "+"), and the B and C layers are actually identical to the A layer only shifted by 
translations of (2a/3, a/3) and (a/3, 2a/3), respectively. Here a is the lattice constant of its unit cell (shown by the black dashed lines). The small (green) and 
large (blue) spheres represent the boron and metal atoms, respectively. 



MoB 4 should be the 2H-type MoB 3 . Meanwhile, the 3R-type MoB 3 
was theoretically presented 32 . It is therefore timely to elucidate the 
long-standing common mystery of this family of TMB 3 and, in par- 
ticular, to provide the essential guidance in the design of superhard 
materials. 

Results 

Polytypism in TMB 3 . The 2H (3R) structure adopts hexagonal 
(rhombohedral) symmetry with space group P6 3 /mmc (R-3m), in 
which boron and metal atoms locate at the Wyckoff 12/ (18/) and 
2b, 2c (6c) sites, respectively. As illustrated in Fig. 1, the boron atoms 
form planar graphitelike sheets (H). The metal atoms sit directly 
above the boron hexagons, skipping one column in every three, 
thus forming honeycomb metal layers (A). The B and C layers are 
actually identical to the A layer only displaced by one atom and two 
atoms, respectively. Accordingly, the stacking sequence AHBH 
(AHBHCH) well describes the 2H (3R) architecture, and both 
consist of identical units of substructure (AH), piled one on top of 
the other in different numbers and in different stacking sequences 
within the unit cells. 

Superhard TMBs have long been believed to exist only in simple 
crystal structures, similar to 2H, 3R or other structure types, follow- 
ing the general idea in the design of superhard materials 2 " 4 . However, 
on the basis of the structure stacking principle mentioned above, we 
may reasonably infer that a very large number of polytypic structural 
modifications of TMB 3 can be built up using alternating metal and 
boron layers by systematically specifying the stacking sequence of 
metal layers. In these architectures, the stacking sequences of boron 
layers are completely the same. The differences amongst various 
structures, however, lie in the stacking sequences of metal layers. 
For the case with two metal layers in the unit cell, the only two 
stacking sequences ( AB and AC) yield an equivalent hexagonal struc- 
ture that is none other than the experimentally observed 2H form 25 . 
Likewise, three metal layers comprise two orders (ABC and ACB) 
that produce a same rhombohedral framework. This framework is 
exactly the theoretically uncovered 3R structure 28,32 . For four metal 
layers, the four stacking sequences (ABAC, ABCB, ACAB, and 
ACBC) equivalently create a previously unreported hexagonal con- 
figuration (called 4H). A complete list of previously unidentified 
structures (called 5H, 6H l5 6H 2 , 7H l5 7H 2 , 7H 3 , 8H l5 8H 2 , 8H 3 , 
8H 4 , 8H 5 , and 8H 6 ) comprising 5-8 metal layers in their unit cells 
can be constructed. In general, the TMB 3 crystals could be imagined 
to display one-dimensional disordered stacking of metal layers along 
the [001] direction. Therefore, polytypism is revealed to be the extra 
hidden degree of freedom in the structure design of superhard WB 3 
and MoB 3 with a large number of polytypes due to the different 
stacking of metal layers amongst the interlocking boron layers. 



Stability of TMB 3 polytypes. The thermodynamic stability of all 
those TMB 3 polytypes can be firmly proved for their viability from 
first-principles calculations (see Supplementary Information for 
details). Based on that, our work is able to systematically identify a 
set of new configurations that are energetically more favorable than 
the experimentally observed but metastable 2H structure. The 
calculated formation energies of fifteen TMB 3 polytypes are 
summarized in Fig. 2. In addition to the concrete values, the 
relative stability of different polytypes for WB 3 and MoB 3 has 
main features in common. The formation energy of the 2H 
structure is found to be the highest amongst the TMB 3 polytypes 
considered, though it is still negative (—1.185 eV/f.u. for WB 3 or 
— 1.256 eV/f.u. for MoB 3 ). This makes the 2H structure metastable 
but thermodynamically viable, which supports the experimental 
observation 25 . Amongst these polytypes, the 3R type has the lowest 
formation energy ( — 1.222 eV/f.u. for WB 3 or —1.283 eV/f.u. for 
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Figure 2 | Formation energy versus the number of metal layers within the 
unit cell for (a) WB 3 and (b) MoB 3 polytypes. The formation energy, 
defined as AE = £(TMB 3 ) - E(TM) - 3E(B), is obtained from the first- 
principles calculations, referring to W and Mo in body-centered cubic 
phases and to B in a-B phase as their respective ground states. The red and 
black dashed lines indicate the lowest and highest energy values, with an 
energy distribution width as small as 37 meV (27 meV) per WB 3 (MoB 3 ) 
formula unit (f.u.), amongst all the calculated polytypes. 
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Figure 3 | Calculated bulk modulus B> shear modulus G, and Vickers hardness inversus the number of metal layers within the unit cell for WB 3 and 
MoB 3 polytypes. In (a)-(c) and (d)-(f), the same symbols denote the same WB 3 and MoB 3 polytypes, respectively. 



M0B3) and becomes the ground state of TMB 3 . Hence, the 2H and 3R 
structures are the extremes in formation energy governing the 
thermodynamic stability of these polytypes. It is interesting to note 
that the energy differences amongst polytypes are extremely small, 
and the maximum separation is as low as 0.037 eV/f.u. for WB 3 and 
0.027 eV/f.u. for MoB 3 . Considering the temperature effect on 
formation energy, such a small energy difference manifests 
polytypism in TMB 3 with a large class of energetically degenerate 
structures. 

As a matter of fact, the newly identified TMB 3 structures can be 
envisaged as different mixtures of the 2H and 3R ones. In these 
layered crystals, the metal layers can be considered as being held 
together by the boron layers. In contrast to the classic layered mate- 
rials (e.g., graphite, /z-BN) that are characterized by weak interlayer 
interactions (e.g., Van der Waals), TMB 3 are dominated by strong 
intralayer B-B and interlayer TM-B covalent bonds 22 ' 28 . Stacking 
faults between the adjacent metal layers, say between AHB and 
AHC, basically do not alter the local environments of B-B and 
TM-B interactions and the energy difference amongst different 
stackings is very small. Hence, stacking faults, resulting from the 
intrinsic short-range behavior of strong covalent bonds, play a 
decisive role in the occurrence of polytypism in TMB 3 . 

The dynamic lattice stability of TMB 3 polytypes can also be sup- 
ported through their phonon dispersion. For the 2H and 3R struc- 
tures, the calculated phonon dispersion curves show no imaginary 
frequencies throughout the whole Brillouin zones, confirming the 
dynamic stability of both basic structures. Since other polytypes are 
viewed as different mixtures of 2H and 3R, the TMB 3 polytypes are 
deemed dynamically stable (see Supplementary Information for 
details). 

Based on the published XRD patterns 11,25 of the experimental 
samples, only the 2H structure can agree with them, since others 
polytypes cannot avoid introducing apparent additional peaks gen- 
erated by tungsten layer which obviously disagrees with the experi- 
mental patterns. Although the newly identified structures are 
metastable with respect to the ground-state 3R phase, they are ener- 
getically more favorable than the experimentally synthesized 2H 
structure and are dynamically stable. Moreover, as mentioned above, 
the energy differences amongst various polytypes are extremely 
small. Considering the temperature effect, the newly identified 



TMB 3 polytypes should be viable. The temperature-dependent 
Gibbs free energies derived by Cheng et at. 25 has revealed that the 
2H phase becomes more stable than the 3R phase above 659 K within 
the GGA method (above 678 K within the LDA method), which 
explains why the ground-state 3R phase has never been observed. 
Therefore, these results provide an extremely useful indication that 
the polytypic TMB 3 structures should be synthesized under appro- 
priate conditions (e.g., at proper temperatures). 

Mechanical properties of TMB 3 polytypes. We also evaluate the 
influence of different polytypic structures on mechanical proper- 
ties of TMB 3 . The obtained bulk modulus, shear modulus and 
Vickers hardness for fifteen TMB 3 structures are displayed in 
Fig. 3. The overall trends are very similar, even though mechanical 
properties of MoB 3 are slightly inferior to those of WB 3 . 
Interestingly, their mechanical properties are hardly affected by the 
number of metal layers in their unit cells, and the differences 
amongst various polytypic structures are also very small. In 
particular, the estimated hardness values (38.0-38.9 GPa for WB 3 
or 36.1-37 GPa for MoB 3 ) do not change much from the calculated 
data with the range of only 1 GPa, regardless of the polytypes 
stacking sequence. The results demonstrated that the intrinsic 
hardness of the TMB 3 polytypes mainly depends on the interaction 
between metal and boron layers. In these layered TMB 3 structures, 
the interlayer TM-B bonds are relatively weaker than the intralayer 
B-B bonds, although they are also strong bonds. This allows the 
layers to cleave readily under a large shear stress. It was verified by 
a theoretical study of the 2H structure that has the lowest indentation 
shear strength along the [110] direction under the (001) plane 27 . 
Therefore, the sliding strength between layers largely limits the 
intrinsic hardness enhancement of the TMB 3 polytypes. 

To further enhance the hardness, the relatively easy sliding of the 
layers under the large shear stress must be suppressed. As mentioned 
above, TMB 3 can form a variety of polytypic structures. Moreover, 
these structures have the nearly same energies. We thus propose that 
one can create a multiphase solid-solution compound with the coex- 
istence of various polytypes. This particular material may include a 
large number of interfaces amongst different polytypic structures 
with different easy sliding directions, and these interfaces would 
strongly hinder the interlayer sliding movement of each polytype, 
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Figure 4 | (a) Low-frequency phonon dispersion curves along the [001] direction of the 2H and 4H polytypes for WB 3 and (b) Estimated room 
temperature lattice thermal conductivity along the [001] direction for WB 3 and MoB 3 . In (a), the changes of Debye temperatures for acoustic modes are 
labeled. 



accordingly allowing the possibility of an extrinsic hardness 
enhancement. A recent experiment 15 showed that a greatly enhanced 
hardness was achieved when pure metals (e.g., Ta, Mn, Cr) were 
added to WB 3 . In contrast to this traditional solid-solution of differ- 
ent compositional substances, the hardening way by creating a mul- 
tiphase compound with various polytypes to restrict the propagation 
of dislocations does not change the chemical compositions of mate- 
rials. Moreover, this mechanism is expected to be applicable to other 
systems. As a matter of fact, the ultrahigh strength (hardness) of the 
recently synthesized nano-twin Cu (c-BN) samples 33 ' 34 has already 
corroborated this. The nature of strengthening is the effective block- 
age of dislocation motion by numerous coherent boundaries, 
although the practical approaches of creating grain boundaries are 
different. 

Thermal conductivity of TMB 3 polytypes. Such a class of special 
superhard materials own ultra-low thermal transport capability, 
completely opposite to the conventional knowledge that 
intrinsically strong chemical bonds in superhard materials would 
lead to high lattice thermal conductivity due to the inherently high 
bulk moduli and sound velocities. The polytypism in TMB 3 greatly 
influences lattice thermal conductivity k l along the [001] direction. 
As well understood, the acoustic part of phonons carries most of the 
heat in a material. The out-of-phase vibrations from a large number 
of randomly distributed atoms within the polytypism- enlarged unit 
cell, together with the inherently polytypism-induced interfaces, 
strongly scatter heat-carrying phonons. Meanwhile, as the number 
of metal layers in the unit cell increases, the Brillouin zone keeps 
"folding-in", and the Debye temperatures 0 D and Debye frequencies 
co D of acoustic phonon modes decrease correspondingly 35 . 
Figure 4(a) plots the phonon dispersions along the [001] direction 
for the 2H and 4H structures of WB 3 . Since the Brillouin zone of 4H 
in this direction becomes smaller, the acoustic modes "hit" the 
Brillouin zone boundary much earlier than the phonons in 2H. 
Hence, w D s are greatly reduced from 4.4 THz and 3.8 THz for 2H 
to 2.7 THz and 2.0 THz for 4H for the longitudinal acoustic mode 
(LA) and transverse acoustic mode (TA), respectively. On the other 
hand, the sound velocities (5180 m/s and 6900 m/s) and Griineisen 
parameters (1.37 and 1.74) for the TA and LA modes remain almost 
the same amongst the 2H, 3R, 4H and 5H polytypes. Given these 
parameters, the lattice thermal conductivity k l could be estimated via 
the Debye- Callaway method 36 . The room temperature k l along the 
[001] direction with respect to the number of atoms in the unit cell is 
plotted in Fig. 4(b). It can be observed that k l decreases with 
increasing number of the random stacking layers, primarily due to 
the reduced Debye frequencies or Debye temperatures. Very low k l 
close to minimal might be achieved after several repetitions of units. 



Such unusual ultra-low lattice thermal conductivity of TMB 3 is very 
unique and could be extremely useful in particular applications, such 
as one-dimensional thermal barrier with superhardness. 

In summary, we show that polytypism exists in superhard WB 3 
and MoB 3 with a large number of energetically degenerate structures 
due to the different stacking of metal layers amongst the interlocking 
boron layers. Although all the polytypic structures have nearly 
equivalently intrinsic hardness, such polytypism may create a mul- 
tiphase solid-solution compound with a large number of interfaces 
amongst different polytypes, and these interfaces will strongly sup- 
press the interlayer sliding movement of each polytype, accordingly 
allowing the possibility of an extrinsic hardness enhancement. 
Furthermore, the polytypic TMB 3 manifest extremely low lattice 
thermal conductivity due to structural disorders and phonon folding. 
Our findings shed light on the structure design of superhard materi- 
als and, at the same time, provide a new strategy in the search for 
superhard materials with additional functionalities. 

Methods 

Calculations on the basis of the density functional theory were carried out with the 
Vienna Ab initio Simulation Package (VASP) code 37 and the electron-ion interactions 
are described by means of projector augmented wave (PAW) with 5d 4 6s 2 (4cf5s 2 ) and 
2s 2 2p l electrons as valence for W (Mo) and B, respectively. The exchange -correlation 
functional within the generalized gradient approximation (GGA) 38 were employed. 
Because of the extremely small energy difference amongst different structures, a large 
cutoff energy of 500 eV and dense Monkhorst-Pack A;-point meshes with a grid of 

0. 02 A -1 for Brillouin zone samplings were adopted to ensure the numerical con- 
vergence of formation energy differences to typically 1 meV/atom. 

Having fully optimized fifteen specific polytypic structures for TMB 3 , their 
formation energies (AE) at temperature T=0K were calculated. Bulk modulus (B) 
and shear modulus (G) are determined by the efficient strain energy method 39 while 
the intrinsic Vickers' hardness is estimated by Chen's model [H = 2(/c 2 G)° 585 — 3, k = 
GAB] 40 ' 41 . Also, the phonon dispersion curves and lattice thermal conductivities along 
the [001] direction for the 2H, 3R, 4H and 5H structures in WB 3 and MoB 3 were 
calculated, (see Supplementary Information for details). 
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